arXiv:1105.3399v1_ [math.GM] 15 May 2011 


An Integral Representation of Kekulé 


Numbers, and Double Integrals Related 
to Smarandache Sequences 


John M. Campbell 
Toronto, ON 
jmaxwellcampbell@gmail.com 


Abstract 


We present an integral representation of Kekulé numbers for P2(n) benzenoids. 


Related integrals of the form [”_ Set) dix are evaluated. Conjectures relating double 


integrals of the form fj” [”, sO dadk to Smarandache sequences are presented. 


1 Introduction 


It is well-known that some characteristics of Kekulé structures may be enumerated in terms 
of recursive sequences [4] [6] [7] [8] [15]. A variety of discrete formulas for A’ (B) (the number 
of Kekulé structures in a benzenoid system B) are known [4]. 

The existence of Kekulé structures in benzenoid systems is commonly considered to be an 
algebraic problem. A well-known discrete identity follows: letting A denote the adjacency 
matrix of a benzenoid system B with n vertices, det A = (—1)?K(B)? (so a benzenoid 
system has Kekulé structures iff the matrix A is non-singular) [5] [8], a result discovered by 
Dewar and Longuet-Higgins (1952) and rigorously proven in 1980 by Cvetkovié (et al.) [4]. 

Inspired by Fourier analysis, we prove an integral equation for Kekulé numbers for pyrenes 
on a string with hexagons removed from both sides (A (P:(n))), and prove an explicit closed 


form evaluation of [*, Sosa dx. 


We conjecture how the related expressions {j” J”, ser) dadk may be used to represent 


(quadratic) continued fractions, and conjecture that double integrals of this form correspond 
to sequences such as the square root of smallest square of the type k?(k+1) and the number of 
solutions to x? =0 mod n, establishing unexpected number theoretic implications. Many of 
the sequences discussed below may be described as Smarandache sequences, such as core(7). 


2 Main Problems and Results 


Given the tetravalency of carbon, it becomes evident that the number of Kekulé structures 
of a benzenoid hydrocarbon is equal to the number of linear factors of the benzenoid system. 
Given the work of Ohkami and Hosoya, the following recurrence formula for the Kekulé 


numbers for pyrenes on a string (i.e. a system of linearly annelated pyrenes, pyrenes being 
a characteristic subunit of all-benzenoid systems) of the form K(P2(n)) holds [4] [15]: 


K(Pr) _ 10.K (Pri) = K (Pas) (1) 


Consider the companion recursion sequence of (1) given by d, = 10dn_1 — dn_2; do = 1, 
d, = 5 [17]. Surprisingly, these recursion sequences may be used to evaluate the following 
integral, letting c, = K(P,): 


™ cos(nz) (—1)"4+1 [2 
—_— dx = ——7(-¢2 dj nj4/= 
x 2+ sin? x . 2 a Cle) + LJ 3) 


Lemmma 1. 


Proof. First consider the non-zero values of the above equation; use induction. Let P(n) 
be the statement “7, sos Ore) dx = (—C, + dny/2)n”. Suppose P(n) holds. So we want 
P(n—2) A P(n—1) => P(n). Verify that P(n) holds for all n < m for some 2 < m. 
We have KP) = 10K (P,-1) — K(P,_2) (1) and dn, = 10d,,_1 — dn_2; do = 1, dy = 5 hab 
So by the assumption of P(n) we have: 


[7 Mes) ae = r0(e5-1 + deaf 


—1 2+ sin? x 


[ COs 28) 9. 2 (—Cno +d [Pe 
n—2 n—2 3 


-- 2+sin?z 
So it remains to prove that: 
L 10 cos(2(n — 1)x) — cos(2(n — 2)x) — cos(2nz) 
—n 2+sin? xz 


It is elementary to evalute the above integral as follows: 


dz =U (2) 


i 10 cos(2(n — 1)x) — cos(2(n — 2)x) — cos(2nz) i 4 sin(2n7) 
—1 2+sin? x a | 


So (2) holds, letting n be an integer (greater than 1). 


Tv cos 


For odd values of n for [”, so da, interpret the integrand as a shifted skew-symmetric 


function. So it remains to prove that ee = ere which is elementary. 
2 2 


Theorem 1. 
m cos(2nz) 


he Dan a 2 
FOP es a + (—1)" cos(2n sin™! V3) 3 


Proof. An alternative evaluation of d,,, from Jasinski (2008) [17] , follows: 


d,, = (—1)" cos(2n sin! V3) 
So given Lemma 1 , Theorem 1 holds. 
O 


Lemmma 2. Let y, be defined by the following integer sequence: Yn = (4k + 2)Yn—1 — Yn-2; 
yo = 0, yr = 4 for some fired k. Let z, be defined by the following integer sequence: 
Zn = (4K + 2)2n-1 — 2n-2, 2 = 1, 21 = 2K +1 for some fixed k. Then: 


™ cos(nx) (—1)"4+1 2 
——_~-d, = oo — nu na —— 3 
ax k + sin? x v 2 ae Ye) + ALB 7) oo 


Proof. First consider the non-zero values of the above equation; use induction. Let P(n, k) 
be the statement, for fixed k, [”, sCmdx = = 1(—Yn + tna} te): We want: P(n — 2,k) A 
P(n—1,k) > P(n,k). So by the assumption of P(n,k) we have: 


[On OD) te = (4k + 2)a(—y a 3 
— 2 Be DT 


k+sin*2z 


(4k+2) cos(2(n—1)2)—cos(2(n—2)x) 
k+sin* x 


a (4k + 2) cos(2(n — 1)x) — cos(2(n — 2)x) — cos(2nz) ie A sin(2n7) 
—n k+sin? x as 
By induction we must also prove that Vk € R(P(0,k) A P(1,k)). 


It is elementary to establish that: 
Q/1+in 
-rk+sin* x 1+k 
So P(0,k) holds for all k. As for P(1,k), it is elementary to establish that: 


So we must prove that [”, cos(2n2) 7 — (), which is elementary: 


(—4Vko/1 + kw + (1 + 2k)(4 + 2i1n(— = 


™ cos(2z) : TTF 
@gkfen se VEVI+FE 


So we must prove: 


) +éIn(z))) 


(—4V7kvV/1 + ko + (1 + 2k) (a + 2i In(— 
VEVI +k 


ia 4 


= (-44 (2k +1) kk+1)" 


es he eT he ee Hin as in( - )) = ane de + (4b 4.20 


& (14 2k)(m + 2iIn(— )+iln(y)) = (4k +2)r & 2i In(—-—==) +iln(; a5 =T 
1+k 1+k 
1+k k , 1+k k , 
ra In(-——) + In ni; 7 Z) =-T1S In(-——) + In(; i 7) = —7i & In(-1) = —77 


So the proof by induction is complete. 


Let f(x,n,k) be an elementary periodic function of the form f : Z* x R? > R given by 
the expression en for fixed integers n and fixed k. 
Clearly f(n,k) is symmetric (about the y-axis). Interpret f(x,n,k) as a shifted skew- 


symmetric function. So we want that f(z +4,n,k) =—f(—v%+4,n,k) . So we want: 


cos((2n + 1)(r+$)) _ _cos((2n ae) 


k+ sin? x k+ sin? x 


= cos((2n + 1)x) cos((2n + )5) — sin((2n + 1)a) sin((2n + 15) = 


—(cos((2n + 1)5) cos((2n + 1)a) + sin((2n + 1)5) sin((2n + 1)a)) 


Recall that n is an integer, so cos((2n + 1)4) = 0. So the proof is complete. 
O 


Reinterpret z, as a multivariable function given by z,, = (4k + 2)Zn-1% — 2n—-2.43 
Zok = 1, zn = 2k +1. Using OEIS one notices that: z,,1 = cosh(2n sinh! 1), a = 
cosh(2n sinh~'(V/2)), 2n,3 = cosh(2nsinh~'(/3)), zn,4 = cosh(2nsinh~'(2)), ... The pattern 
is clear. 

Reinterpret y,,, as a multivariable function given by yn,4 = (44+2)Yyn—1,4—Yn—-2,k) Yo,k = 9, 
Yi, = 4. Using OEIS one notices that: 


((3 + 20/2)" — (3 — 20/2)") 


Yn = 


((5 + 2/6)" — (5 — 2/6)") 


Given the above results one deduces that: 


Croat =(7=941/719)) 


Yn2 = 


Yn3 = SS 
5 


Also one may deduce that: 


iL n n 
tna = pel((9 + 2V2H)" — (9 — 2V20)") 


So we conjecture: 


Lemmma 3. 


1 
= — (ORD ee DSR S1= 2 he 1)” 
Unk = va(k +1) Vek +1) 


k(k+1 


Zn,k = cosh(2n sinh” '(vk)) 


Proof. We have Yn, = (4k +2)Yn—1,b — Yn—2,k) Yo,k = 0, Yre = 4, and 2n4 = (4k + 2)Zn 1p — 
2-08 Me = 1 tig Se 1 


Let P(n,k) represent the statement “ yn~ = Voorn (2% +1+2,/k(k + 1))"-(2k+1- 


2\/k(k +1))")”, for some fixed k. Verify that P(0,k) and P(1,k) hold: 


1 
Yo,k = s~——((2K +1 4+ 2/k(k + 1))° — (2k +1 — 24/k(k + 1))”) 
vee +1) 


k(k +1) 
Vik = as +1+2/k(k+1))' — (Qh+1—24/k(k+1))") 


So by the assumption of P(n,k) we want P(n —2,k) A P(n—1,k) = P(n,k), and have: 
1 
(4k + 2)Yn—14 = (4k + a +142 /k(k + 1))" 7? — (2k +1—2)/k(k + 1))*4) 
+ 
1 
Yn—2,k = t———— ((2k +14 2/k(k + 1))"-? — (24 4+1—-2/k(k + 1))"?) 


k(k +1) 


So it remains to prove that: 


(26 41 42/e(k +1)" - (2k +1-2,/e(R + 1)") = 


k(k +1) 


1 
ke + 2) —————((2k +14 24/k(k n-l _ (94 2, Ee n-1 
(4k + 2) Teste +1+ (k + 1)) (2k +1—2\/k(k+1))"~) 
il 
—————— ((2k + 14+ 2\/k(k + 1))"? — (2h + 1 — 24/k(k +1))”?) 
k(k +1) 


(2k +14 2,/k(k+1))” — (2k +1-2/k(k+1))” = 
(4k +2)((2k +14 2/k(e+ 1)" — Qk +1 -2fe(k+1))") 
—((2k+1+2/k(k + 1))"? — (2k +1-2V/k(k + 1))”’) 


Simplify the right-hand side expression and the above equation holds. 
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As for Zn, we have: Zn, = (4k + 2)Zn_1,4 — 2n—2,k3 20k = 1, 21,4 = 2k +1, and we want 
Zo, = land 2, = 2k+1 


zo,h = cosh(2 x Osinh~!(Vk)) 
21," = cosh(2 sinh '(vVk)) 
Expand and simplify, and the above holds. 
(4k + 2)zn—14 = (4k + 2) cosh(2(n — 1) sinh7!(Wk)) 


Zn—2,k = cosh(2(n — 2) sinh7'(Wk)) 


So it remains to prove that: 


cosh(2n sinh~!(Wk)) = (4k + 2) cosh(2(n — 1) sinh~!(Wk)) — cosh(2(n — 2) sinh!(Wk)) = 
(4k + 2) cosh(2nsinh7'(Vk) — 2sinh!(Vk)) — cosh(2n sinh (Wk) — 4 sinh!(Vk)) 


Expand and simplify, and the proof holds. 


O 
Theorem 2. 
™ cos(nz) —_ 
_ck+sin?z 
—1)" 2k +1+4+2/k(k +1))'2! — (2k4+1-2,/k(k + 1))!2! 
(ay s1 | Gk bay hier Dye — k-+ 1 2y Rie rD) 
2 k(k +1) 
2 
+ cosh(2 [2 | sinh7?(/k)) —————) 
2 k(k +1) 
Proof. The above theorem follows from Lemma 2 and Lemma 3. 
O 


Note that k& is not restricted to integers. Integrating the above with respect to k yields 
strange results, often related to Smarandache sequences. 
Consider the following function for even integers n: 


1 cos( nx) 
5 4 i"\(2 nf - wt 
I(r) = rian a rt+sin?x re 


Consider: 


f(10) = —41(—41 + 29,/2)x 


f(12) = —140(—140 + 99/2) 
f(14) = —239(—239 + 1692) 
f(16) = —816(—816 + 577V/2)x 


f (18) = —1393(—1393 + 9852) 


Remarkably the sequence ..., 41, 140, 239, 816, 1393, ... corresponds to the numerators 
of the lower principal and intermediate convergents to \/2 [9], and the sequence ..., 29, 99, 
169, 577, 985 , ... corresponds to the denominators of the lower principal and intermediate 
convergents to V2 [1]. 


Conjecture 1. Let A, represent the numerator of the nth convergent to V2 and let P, 
represent the nth Pell number, i.e. the denominator of the nth convergent to V2. 
If 5 is even then: 


1 
a(-5 +") +i") af [ dadt = 2P,)9(—2Ppj2 + AnjoV2)m 
If 5 is odd then: 


il cos(nx 
we ae Need) nf a dt = Anpo(— Anja + Pappv2)n 


(Note that the double integral [', {7, cote) dadt displays similar behaviour.) 
Moreover, 


Conjecture 2. For all integers m, and for all even integers n, the expression 


cos(nz) 
[ [ didt 
rt+sin? x 
is equal to a number of the form a(b+ cVdm)m, where dm is equal to the mth smallest k 
such that k?(k +1) is a square, and a, b, and c are rational. 


The aforementioned sequence (2, 6, 3, 5, 30, 42,...) is given in [12]. To illustrate the 
above phenomenon, let n = 4, and evaluate the above double integral for m = 1, 2,3,..., and 
note the appearance of the aforementioned integer sequence (2, 6, 3, 5, 30, 42, ...): 


4(—4 + 3/2) 
4(—12 + 5/6) 
8(-12 + 7/3) 
8(—20 + 9/5) 
4(—60 + 11/30) 
4(—84 + 13/42) 
8(—56 + 15\/14)x 
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Another example of the above double-integral-quadratic-convergent phenomenon: 


Conjecture 3. Let C,, represent the numerator of the nth convergent to V3 and let Dy 
represent the denominator of the nth convergent to V3. Ifn=2 mod 4, then: 


cos( nx) 
zn Inf facet = Ca (-Cy + Dg V3) 


at+sin?x 


(A similar phenomenon holds for n = 0 mod 4.) 
One may easily conclude that the above phenomenon does not hold for all rational m; 
consider: 


Conjecture 4. For all even positive n, 


cos(nz) 4x 32-4 
Nig gO eee 
t+sin? x 32 
Further conjectures related to number theory follow: 


Conjecture 5. Consider the sequence of numbers given by: 


cos(22) 
I a dxdt 
rt+sin2 2 
Factoring the maximum integer from these expressions, the resulting coefficient of the 
square root is given by the sequence of the number of solutions to x2 =0 mod n. 


This sequence (1, 1, 2, 1, 1, 1, 2, 3, 1, 1, 2, 1, 1, 1, 4, 1, ...) is given by [18]. Note 
the appearance of this sequence (and of [12]), evaluating the above double integral for m = 
1,2,3, 20: 

4(—1+ /2)n, 4(-24 V6)a, 4(-3 + 2V3)z, 8(-2 + /5)a, 4(—5 + V30)z, 4(-6 + /42)z, 
A(—7 + 24/14), 8(—4 + 3\/2), 12(—3 + 10), 4(—10 + 110) x, 4(—11 + 2,\/33)z, 8(— 
/39)r, 4(—13 + /182)m, 4(—14 + 210), 4(—15 + 4/15), 16(—4 4+ V17)z, 

Factoring the maximum integer from these expressions, note the sequence of (the absolute 
value of the) initial negative coefficients (1, 2, 3, 2, 5, 6, 7, 4, 3, 10, 11, 6, 13, 14, 15, 4, ...), 
which corresponds to the Smarandache sequence [10]. We conjecture that: 


Conjecture 6. Consider the sequence of numbers given by: 
cos(22) 
I a dicdt 
rt+sin2 x 
Factoring the maximum integer from these expressions, the (absolute value of the) result- 


ing negative coefficient is given by the sequence given by the smallest number whose square 
is divisible by n: (1, 2, 3, 2, 5, 6, 7, 4, 3, 10, 11, 6, 13, 14, 15, 4, 17, 6, ...). 


Conjecture 7. Consider the sequence of numbers given by, form = 1, 2,3,...: 


= = cos(2z) 
xt+sin?x ra 
Factoring the maximum integer from these numbers, the (absolutel value of the) resulting 
negative coefficient is given by the smallest number whose square is divisible by n. 
Again the Smarandache sequence [10] arises. 
Consider the sequence of numbers given by, for m = 1, 2,3,...: 
welt 


=a 1 
[ ent 
et sin? x 


Unexpected patterns arise in . evaluation of the above double integral. Let: 


mit 
2 — 
a [ eee adit = —n(In(a2,) — 2 1n(Bm + Cm/dm)) 


Using Sloane’s online sequence recognition tool, we conjecture that: 


Conjecture 8. b,, corresponds to A165367, the trisection a(n) = A026741 (3 «n+ 2), 
where A026741 is n if n odd, n/2 if n even. a?, corresponds to A168077, which is 


a(2n) =A129194(2n)/2; a(2n+1) = A129194(2n +1), where A129194 is n?(2 — (—1)*). 


We presently leave as open problems the evaluation of c,, and d,, (defined above). 
Let: 


m™m (42) 
[ es om ———drdt = Tam(—bm + Cn\/ dm) 


rt+sin? x 
Using Sloane’s online sequence recognition tool, we conjecture that: 


Conjecture 9. c,, corresponds to A165367, the trisection a(n) = A026741 (3 «n+ 2), 
where A026741 is n ifn odd, n/2 if n even. 


Let: 


mxV/2 pr 1 


Using Sloane’s online sequence recognition tool, we conjecture that: 


Conjecture 10. b,, corresponds to the Smarandache sequence AO07913, which is the 
squarefree part of n: a(n) = smallest positive number m such that n/m is a square. 


Note that the Smarandache sequence corresponding to b,, is also referred to as core(m) 
[11]. We thus have related Fourier-type integrals to squarefree numbers (through [11]), given 
that the above double integrals may be represented as follows, by Fubini’s theorem: 


ia (f" cos(nz) dee)dt 
rt+sin?x 
= a cos(nx)(—In(sin? x) + In(m + sin? x))dx 


Conjecture 10 may easily be extended to all non-perfect-square integers, with the Smaran- 
dache sequence A007913 [11] again arising. 

We presently leave the above conjectures as open problems. A closed form evaluation 
ie be ea ost dcdt would be ideal, especially considering the multifarious connections to 
nimber theory. 


3 Other Results from Fourier Analysis 


Consider elementary 27-periodic functions of the form f : Q? x R > C given by f(x) = 
In(sin(ax)+cos(ba)), where a and bare rational. Letting a = 1 and b = 1 the evaluation of the 
sine and cosine Fourier coefficients of f(2) results in expressions involving 2F\(1, —$,1—4,1), 
which is not defined for positive even integers n. We evaluate the aforementioned coefficients 
as follows: 


Proposition 1. 


Yn € arf cos(nz) In(sin z + cos x)dx = 


=H 


e ae] 1 n +3) mod 4 n 27 
yee oP +t 1 4 eee iain 
Proposition 2. 
” =] | m5 | _1)n+3 Rec (n+5) mod 4 =| | 24? | 
Yn € z(f sin(nz) In(sin x+cos «)dx = ee ee eee 


Note that the former series may be expressed using the Kronecker symbol and related 
modular forms [19], and that | Setsimod* | (—1)!4! corresponds to the inverse of the 8th cy- 
clotomic polynomial [13]. 

Given Proposition 1 and Proposition 2, determine the Fourier series for In(sin x + cos x): 


In(sin x + cos x) AG * ae 
+ = cos(nz)((—1)!" a + er simeds (—1)l4 =) 


1) L222) (pnt 4 1)j «| et mod 4 [242 | 
epi ar 


Van 2B] +1 


Letting x = 7, use the above technique to prove: 
mS EDIE (yt +0) 
V2 n=1 2 LS! or 
n mod n 
| ets) mod t J(—1)t4! 
n 


To conclude, in addition to presenting the above two series, we have rigorously proven the 


new evaluation of {7, en de, and have presented unusual connections between Fourier- 


type double integrals and Smarandache number theory. 


4 Appendix 


Mathematica function for Theorem 2: 


Theorem2 = 
Simplify [FunctionExpand[ (-—(((2#1 + 1 + 2 Sgqrt[#l1 (#1 + 1)])% 
Floor[#2/2] - (2 #1 + 1 - 2 Sqrt[#1 (#1 + 1)])% 


Floor [#2/2])/Sqrt [#1 (#1 + 1)]) + (Cosh[ 
2 Floor[#2/2] ArcSinh[Sqrt[#1]]]*(2/Sqrt [#1 (#1 + 1)])))* 
Pix ((C(=1) *#2) + 1)7(2)) 1] «& 


Integration results for Conjecture 8: 
—n(In(4) — 2In(7 + 3/5)) 


—r(In(2 — 2In(17 + 266), 
—7(In(9) — 2In(10 + ve 


n(4 


5) 

9) 

—m(In(49) 
~m(In(16) 

l 

( 


30)) 
7) 
30 
=) 
—r(In(36) — 2In(19 + 5 /15)) 
—r(In(169) — 21n(41 + 6/42) 
—r(In(49) — 2 1n(22 + /435)) 
—n(In(225) — 2In(47 + 8/31), 
—n(In(64) — 21n(25 + a 61)) 

70) 

3)) 


—r(In(361) — 21n(59 + 4\/195)), 
—r(In(100) — 2In(31 + +/861)), ... 


Lil 


As indicated above, the sequence 4, 25, 9, 49, 16, 81, 25, 121, 36, 169, 49, 225, 64, 289, 
81, 361, 100, ... corresponds to A168077 [2], and the sequence 7, 17, 10, 23, 13, 29, 16, 35, 
19, 41, 22, 47, 25, 53, 28, 59, 31, ... corresponds to A165367 [3]. 


Integration results for Conjecture 9: 


4(-12 + 56)r, 
2(-15 + 4y/15)rr, 
§(—28 + 11V/7)n, 
3(-15 + 7V5)r, 
24 17/66), 
(—91 + 10/91)r, 
(—120 + 23,/30)z, 
(—51 + 13V17)r, 
(—380 + 29,/190)r, 
#(—231 + 16,/231)z, 
= (—276 + 3569)z, 
3 (—65 + 19/13), 
72 (—252 + 41/B)r, 
ig(— 435 + 22,\/435)r, 
se (—248 + 47V/31)r, 
w(— 561 + 25./561)z, 
sa (—420 + 53V/70)rr, 
2(_ 703 + 28./703)7, 
aor (— 780 + 59\/195)r, 
= (—861 + 31/861)z, 


<8 


As indicated above the sequence 5, 4, 11, 7, 17, 10, 23, 13, 29, 16, 35, 19, 41, 22, 47, 25, 
53, 28, 59, 31 corresponds to A165367 [3]. 


Integration results for Conjecture 10: 


Qn In(1 + 2/2 + 21/2 + V2), 
Qn In(1 +472 4+ 21/2(4+ V2)), 
2m In(1 + 6V2 + 2,/3(6 + V2)), 

Qn In(1 + 8/24 4/84 V2), 

Qn In(1 + 10/2 + 2\/5(10 + V2)), 
Qn In(1 + 12./2 + 24/6(12 + /2)), 
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19, 


Qn In(1 + 14/2 + 24/7(14 + V2)), 
Qn In(1 + 16/2 + 41/2(16 + /2)), 
Qn In[1 + 18/2 + 6/18 + V2)), 
Qn In{1 + 20/2 + 24/10(20 + /2)), 
Qn In[1 + 22/2 + 24/11(22 + +/2)), 


Qn In{1 + 26/2 + 2,/13(26 + /2)), 
Qn In[1 + 28/2 + 2,/14(28 + /2)), 
Qn In{1 + 30/2 + 24/15(30 + /2)), 
27 In 
Qn In[1 + 34/2 + 24/17(34 + /2)), 
Qn In[1 + 36/2 + 61/2(36 + V2)) 
Qn In[1 + 38/2 + 24/19(38 + /2)), 
Qn In[l + 40/2 + 44/5(40 + /2)), ... 


wary 
+ 
iS 
a 
+ 
ee) 
& 
+ 
3 


As indicated above, the sequence 1, 2, 3, 1, 5, 6, 7, 2, 1, 10, 11, 3, 13, 14, 15, 1, 17, 2, 


5, ... corresponds to the Smarandache sequence A007913 [11]. 
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